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It is demonstrated that the claim in the paper "Crossing uj = —1 by a single scalar field on a 
Dvali-Gabadadze-Porrati brane" by H Zhang and Z-H Zhu [Phys.Rev.D75,023510(2007)], about a 
prove that there do not exist scaling solutions in a universe with dust in a Dvali-Gabadadze-Porrati 
(DGP) braneworld scenario, is incorrect. 
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Since the discovery that our universe can be cur- 
rently undergoing a stage of accelerated expansion 
many phenomenological models based either on Einstein 
General Relativity (EGR), or using alternatives like the 
higher dimensional brane world theories have been 
invoked (for a recent review on the subject see reference 
3). The latter ones, being phenomenological in nature, 
are inspired by string theory. 

One of the brane models that have received most at- 
tention in recent years is the so called Dvali-Gabadadze- 
Porrati (DGP) brane world 4].^ This model describes 
a brane with 4D world-volume, that is embedded into a 
flat 5D bulk, and allows for infrared/large scale modifi- 
cations of gravitational laws. A distinctive ingredient of 
the model is the induced Einstein-Hilbert action on the 
brane, that is responsible for the recovery of 4D Einstein 
gravity at moderate scales, even if the mechanism of this 
recovery is rather non-trivial Q . The acceleration of the 
expansion at late times is explained here as a consequence 
of the leakage of gravity into the bulk at large (cosmolog- 
ical) scales, so it is just a 5D geometrical effect, unrelated 
to any kind of misterious " dark energy" . 

The study of the dynamics of DGP models is a very 
atractive subject of research. It is due, in part, to the 
very simple geometrical explanation to the " dark energy 
problem" , and, in part, to the fact that it is one of a very 
few possible consistent infrared modifications of gravity 
that might be ever found. In particular, there can be 
found studies of the dynamics of a self-interacting scalar 
field trapped on a DGP brane by invoking the dynami- 
cal systems tools, which have been proved useful to re- 
trieve significant information about the evolution of a 
huge class of cosmological models. In this regard, the ex- 



ponential potential represents a common functional form 
for self-interaction potentials that can be found in higher- 
order [3] or higher-dimensional theories H. These can 
also arise due to non-perturbative effects [9| . 

A dynamical study of DGP models with a self- 
interacting scalar field trapped on the DGP brane has 
been undertaken, for instance, in reference [10| for an ex- 
ponential potential, to show that crossing of the phantom 
barrier w = — 1 is indeed possible in DGP cosmology with 
a single scalar field (see also [lH in this regard). However, 
the authors of that paper do not study in detail the phase 
space of the model and, in correspondence, they are not 
able to find critical points. Their claim that scaling so- 
lutions do not exist in a universe with dust on a DGP 
brane (only the Minkowski cosmological phase is consid- 
ered), seems to be in contradiction with known results. 
In fact, in the 4D limit (the formal limit when, in the 
DGP model, the crossover length Tc = ^1/^^ oo)^ the 
results of reference [l2| have to be recovered, or, at least, 
approached, since the investigation in [l^ is just a gen- 



eralization of the one reported in [12l | , to include higher- 
dimensional behaviour dictated by the DGP dynamics. 
Even if one expects that the late-time structure of the 
phase space is modified by the contribution of the DGP 
brane, these modifications should be associated with the 
stability of the corresponding critical points rather than 
with their mere existence. 

In the present comment we show which is the source 
of the incorrect result of reference p^ . and we perform 
an exhaustive analysis of the phase space for the DGP 
model with a self-interacting scalar field trapped on the 
brane, by using the same variables of [To| . It is shown, in 
particular, that there is actually an isolated critical point 
associated with the matter-scaling solution, even if it is 
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always a saddle point in phase space. 

The starting point is the Friedmann-DGP equation on 
the brane (equation (5) of reference [lol|): 



These variables are subject to the Friedmann constraint 
(equation (27) in [Ifl]) coming from the Friedamnn-DGP 
equation H]): 



p + Po + Opo 



1 + ^ 
Pq 



1/2 



(1) 



where H — a /a is the Hubble parameter, a is 
the scale factor, k is the spatial curvature of the 
three-dimensional (maximally symmetric) Friedamnn- 
Robertson- Walker (FRW) space ~ taken here to be van- 
ishing: k — Q - , and 9 — ±1 denotes the two branches 
of the DGP model (the two possible ways to embedd the 
DGP brane into the Minkowski bulk). In what follows 
we shall consider the Minkowski cosmological phase, i. 
e., the case = — 1, exclusively. The total energy density 
on the brane p, includes dust matter and the scalar field 
"fluid": 



+ f + h'^ - 
or, alternartively 

f = l-: 



1/2 



±V2b. 



1, (8) 



(9) 



In what follows we shall consider only the "-I-" sign in 
(|9|), since we are focused here on expanding universes 
only, while the oposite sign "— " corresponds to contract- 
ing universes. Thanks to the constraint ([9]), the dimen- 
sion of the phase space can be reduced from 4 to 3. The 
consequence is that, only 3 of the 4 ordinary differen- 
tial equations (22-25) in [T3|, are independent. Here we 
choose the following as independent equations: 



P = P<j}+ P-n 



(2) 



The effective "density" po relates the strength of 5- 
dimensional gravity with respect to the 4-dimensional 
gravity. 



(3) 



where, as usual, is the crossover radius. It is evident 
that 4-dimensional behavior is associated with the formal 
limit Pq — > 0. In the model of interest, 



+ p^ = \^'' -v{4>), (4) 

where a dot denotes derivative with respect to the cosmic 
time t and V{(j)) is the self-interaction potential, taken 
here to be in the form of a single exponential: 



V{(l>) ^ Vo e-^f-^/^). 



(5) 



Here A is a constant parameter and Vq denotes the initial 
value of the potential. 

In order to write the equations of the present model: 
the Friedmann equation ([T]) plus the continuity equations 



X 



^ay{2x' + l'') + ^\xy, 



y = 

h' = -^ah{2x'+F), 



where 



1 



1 



62 



-1/2 



(10) 

(11) 
(12) 

(13) 



and the prime denotes derivative with respect to the time 
variable s = In a. As already said, due to the constraint 
([8|), the variable I can be written as a function of the 
remaining variables x, y, and b (see equation so 
that the differential equations (|10m2[) can be written in 
the alternative form: 



-ax{l + - + V2b) + 3, 



x^^Xy^ilA) 



y' = -'^ay{l + x^~y^ + V2b) + ^Xxy, 



b' = --ab{\^x^ 



y 



V2b), 



(15) 
(16) 



where, now 



Prn + 3Hp„i = 0, p^+ iH{p^ + p^) = 0, 



(6) 



in the form of a (autonomous) dynamical system, the 
following dynamical variables are chosen (see equations 
(18)-(21) in [H): 



a = 1 - 



Vb^ + 2V2b + 2 



(17) 



The first step is to identify the phase space for our 
model, which is given by the non-compact 3-dimensional 
region:'^ 



V 



V6pH' VS/ii?' VS^H' VipH' 



(7) 



Our focus will be on expanding universes only. 



3 



TABLE I: Critical points of the autonomous system of differential equations (|10ll2p . 



Pi X y b Existence ujrf, q 



Pi 











All A 


undefined 


1/2 


1 




±1 








All A 


1 


2 





P-i 


[31 

V 2 A 


/i" 1 
V 2 A 





A^ >3 





1/2 


1 " A^ 


P4 


A 







A^ <6 


^-1 

3 


^-1 






* = {{x, 5) : < + < 1 + %/2&, < y, < 6}, 

(18) 

Note that, a straightforward analysis of the ordinary 
differential equations in the above autonomous system of 
equations pom2p . shows that we can not have x = y — 
b = I = simultaneously, since this would imply that the 
constraint ^ is not obeyed. Therefore, the critical point 
found in the reference ^IQ] , where x = y = b = l = 
at the same time, does not really belong in the phase 
space of the model under study. By the same reason, 
the other point found in [lo| : x — y ~ I ~ 0,h ^ const, 
does not belong in ^ neither. In fact, from the constraint 
^ it follows that, if x = y = Z = 0, then h = -1/V2, 
which is not in 5* since we are considering expanding 
FRW universes only. 

Going a step forward we can realize that the only crit- 
ical points of the autonomous system p4j|16p are associ- 
ated with the four-dimensional limit ^ 00 /Oq — > 
6 = 0.^ This case coincides with the one studied 
in reference [l^]. The constraint ^ translates now into 
the following relationship: 

f ^l-x^ ~y^, ^ a = l. (19) 
We are led with the two-dimensional system of equations: 

x'^~\x{l^x^-y^)+ix-^\y\ 

y' = -\y{\ + x^ -y^) + ^\xy. (20) 

The critical points of (|10m2|) are summarized in the ta- 
ble ID These coincide with the ones found in reference 
as it should be. The point Pi = {x,y,b) = (0,0,0) cor- 
responds to the matter-dominated solution (recall that. 



* The critical points {x,y,b) = (0,0,-1/^2) and (0,0,-0.91 ± 
i 0.68) do not belong in the phase space ^ so that these are not 
critical points of the present model. 




FIG. 1: Trajectories in phase space for different sets of ini- 
tial conditions. Pictures in the upper part of the figure 
show different levels of magnification. The fiux in time t is 
shown in the left-hand picture in the lower part of the figure. 
The value of the parameter of the exponent in the potential 
V — Voe^'^^^, has been chosen to be: A = 2. In the right-hand 
picture in the lower part of the figure we show the projection 
onto the plane {x,y,0). The liinetic energy dominated so- 
lution P^ is the past attractor. The scalar field-dominated 
solution seems to be the future (late-time) attractor. This 
is just a "mirage" due to the fact that the brane effects are 
ignored in the projection. 



in this case P — 1): = pm/5ii^. The kinetic energy- 
dominated solutions (x^ = I, p = ^ H"^ = cjji'^/Qp'^) 
correspond to the points in TabU These are always 
unstable source points in phase space, i. e., the past 
attractor for every trajectory in the phase space. 

The remaining critical points are the scalar field dom- 
inated solution (point P4) and the matter-scaling solu- 
tion P3. In the first case, since P — — 0, then 
x^ +y^ = 1 =^ ij2 = (0V2 + V)/3tP. This solu- 
tion exists whenever A^ < 6 and is accelerating if A^ < 2. 
The matter-scaling solution (critical point P3) exists for 
A^ > 3. In this case x = y, P = 1 — 2x^. The scalar field 
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FIG. 2: Trajectories in phase space for different sets of initial 
conditions, for A = 3 in the exponent (upper part of the fig- 
ure). In the lower part of the figure we show the flux in time 
(picture to the left), and projection on to the plane {x,y,0) 
(picture to the right). Due to the ignorance of the brane ef- 
fects in the projection, the matter-scaling solution seems to 
be the late-time attractor. 

fluid mimics dust matter (w^ = = OY 

Unlike the standard result in both solu- 

tions (scalar field-dominated and matter-scaling critical 
points) represent saddle points in the phase space of the 
DGP model. Looking at the projections of the phase 
space onto the plane (x, y, 0) - the plane associated with 
four-dimensional behavior in the lower figures in Fig. [T] 
and Fig. [2]- so that we are ignoring higher-dimensional 
effects, one can think that these solutions, depending on 
the values of the free parameters, can be late-time attrac- 
tor points in phase space (the standard result in [13] )■ 
However, since in the DGP model higher-dimensional ef- 
fects modify the late-time dynamics, the actual situa- 
tion is very different. At late times the trajectories in 
phase space leave the plane {x, y, 0) and asymptotically 
approach to increasingly large values of the parameter 



b - associated with the DGP brane effects - and of the 
scalar field kinetic energy so that there is no iso- 

lated critical point in the phase space that could be as- 
sociated with a unique future (late-time) attractor. The 
phase space pictures in the figures [T] and [2l show this very 
interesting fact. It is apparent that the phase space tra- 
jectories leave the (x, y, 0)-plane at different places (and 
times), and these approach to different regions of the 3D 
phase space. It is clear, also, that the critical points as- 
sociated with the matter-scaling solution, and with the 
scalar field-dominated solution, can be only saddle points 
in the phase space of the DGP cosmological model. 

We can conclude this comment by identifying the 
source of the incorrect result reported in reference [lo| . 
regarding the existence of matter-scaling critical points 
in a dust universe in the Minkowski cosmological phase 
of the DGP model with a scalar field trapped on the 
brane: the inaccurate identification of the phase space 
corresponding to the model. The existence of matter- 
scaling solutions is expected from the beginnig since, in 
the 4-dimensional limit when standard Friedmann behav- 
ior is recovered, we are left with the case studied in the 
reference [l^, where these solutions were identified as 
critical points in phase space. That is true even if it is 
expected that the stability of these points is modified by 
the infra-red (DGP) brane effects. The kinetic energy- 
dominated solution is always the past attractor (as in 
[l^ ) since, at early times, the brane effects can be safely 
ignored so that the standard cosmological dynamics is 
not modified. 

A graphic illustration of the above discussed features is 
given in the figures[l]and[2]above, where the phase space 
trajectories uncover the main features of the dinamical 
system of interest in [10]. 
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